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O . Abstract 

m 

The relation between the heavy-quark field in QCD and the corresponding field 
' r~| ', in HQET is derived up to three loops, and to all orders in the large- /3o limit. The 

Qh| corresponding relation between the QED electron field and the Bloch-Nordsieck one 

Q^' is gauge invariant to all orders. We also prove that the MS anomalous dimension of 

^ . the QED electron field depends on the gauge parameter only at one loop. 

in 

^ ■ QCD problems with a single heavy quark Q can be treated in a simpler effec- 

|0 . tive theory — HQET, if there exists a 4-velocity v such that the heavy-quark 

momentum is p = mv + k (m is the on-shell mass) and the characteristic 
residual momentum is small: k <^m. QCD operators can be written as series 



(N 



O I in 1/m via HQET operators; the coefficients in these series are determined by 

^ ■ matching on-shell matrix elements in both theories. For example, the heavy- 

light quark currents have been considered at the leading (zeroth) order in 1/m 
• r-j , up to three-loop accuracy [II2]- 

rS 

' Here we shall consider the heavy-quark field Q. Though its matrix elements are 

not directly observable, its matching coefficient can have some applications. 
For example, it is not possible to simulate heavy quarks on a lattice directly, 
because at present we cannot have lattice spacings a <C 1/m. On the other 
hand, simulating HQET on a lattice only requires a ^ 1/Af^. It is possible 
to obtain the HQET heavy-quark propagator in the Landau gauge from such 
simulations. Then, if we know the matching coefficient, we can reconstruct a 
fundamental QCD quantity — the heavy-quark propagator as a function of x. 

At the tree level, Q is related to the corresponding HQET field Qy (satisfying 

^Qv = Qv) by mm 
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The matrix elements of the bare fields between the on-shell quark with mo- 
mentum p = mv + k and the vacuum in both theories are given by the on-shell 
wave-function renormalization constants: 



<0|Qo|Q(p)> 



/70s 



1/2 



u{p), <0\Qm\Q{p)> 



(HQET renormalization constants are denoted by Z). The Dirac spinors are 
related by the Foldy-Wouthuysen transformation 



u{mv + k) 



2m 



Uv{k) 



Therefore, the bare fields are related by 



V 2m / Vm^ 



(3) 



where the bare matching coefficient is 



(4) 



(we use the covariant gauge: the gauge-fixing term in the Lagrangian 
is —{d^AQ^)/{2aQ), and the free gluon propagator is {—i/p'^){g^iu — (1 — 
ao)p^Pu/p'^)', the number of fiavours in QCD is Uf = rii + 1). The C(l/m) 
matching coefficient in (|3]) is equal to the leading one, zq; this refiexes the 
reparametrization invariance [7]. The MS renormalized fields are related by 
the formula similar to ([3]), with the renormalized decoupling coefficient 



ZQ(a("')(/i),a("')(/i)) 
ZQ(a^-^^)(/i),a("^+i)(/i)) 



(5) 



If there are no massive fiavours except Q, then Zq = 1 because all loop cor- 
rections are scale-free. The QCD on-shell renormalization constant Zq con- 
tains the single scale m in this case; it has been calculated [8] up to three 
loops. The three-loop MS anomalous dimensions of pi9] and Q p^fTT] 
are also known. We have to express all three quantities ZQ{gl^'^^\a'^^~^^^), 
ZQ(a("'+i)(/i),a("'+^)(/x)), ZQ(a("')(/i),a("')(/i)) via the same variables, say, 
a^'^'\ij,), see p2]. An explicit formula expressing a^J^'-~^^\fx) via a^"'^(/i) 
and L = 21og(/i/m) (m is the on-shell mass) can be found in [13]. The corre- 



2 



spending relation between a("'+^)(//) and a^^'^n) is 



1 



3^ + ^^^ + 9 + 



T, 



47r 



13 



7r2-13 



Stt^ + 169 
12 



47r 



+ 



(6) 



Our main result is the renormalized matching coefRcient 

z{i^) =1 - (3L + 4)Cf^^^ + {Z22L'' + Z21L + Z20) Cf i^^^'^^^^^ 



47r 



+ (z^^L^ + ^32^2 + -ZaiL + ^3o) Cf ^ 



47r 



+ 



where 



^22 



Z2\ 



Z20 



Z32 



9 11 
^-Cf-y'^a + ^TfUi, 

'^^r 215 „ ^38 

/ 241 \ 

= ( IGtt^ log 2 - 24C3 - IStt^ + — j Cf 

+ (-Stt^ log 2 + I2C3 + Svr^ - ™) Ca 

/4 2 113\^ / 16 2 947\ 
V3 6y 'V3 18/ ' 

= - ^C^F + ^CfCa - ^Cl - QCfTfUi 



■Cl - 42CFTFni - ACfTf 



45 

= - yC| + 135CpCa 



2^31 = 



7831 

-487r^ log 2 + 72C3 + 397r^ - — d 

8 



+ 



"424 



-r 



331 



^log2 - 224C3-^7r2 + 



13307' 
24 



(7) 



3 



+ 



2^9^ 1 
"45" +4^3 + 4l« 



(ni) 



176 2, 325^ 2 4 110 2 73981 



^2 



+ 



128 2, ^ 92 2 613\ 

— TT^ log 2 + I6C3 + ^TT^ - — 1 
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CfTfUi + IGtt^ - 



1013\ 



/64 2, . 32 2 10816\^ ^ 

+ ( y TT^ log 2 + I6C3 - -g-TT^ + \ CaTfUi 



(121 , 352 2 11278\^ ^ 



^30 



32 2 1336 
224 



128 
~9" 



-TT 



3908^ 
3568 



20. 
"9 



-1792a4 - — ^ log^ 2 + 967r^ log^ 2 + 

o 2. -,c..n. 76 4 4801 2 
+ STT^Ca - 1256C3 - — TT^ - -^tt' - 
15 9 



TT^ log 2 - 2OC5 

3023] . 
12 ^ 



+ 



32 4, 4 1448 2, 9 2752 
--04 - - log^ 2 - -^TT^ log2 2 - -^TT^ log 2 + 580C5 



1807rX3 



2312 



-C3 + 



6697 



+ 



7 A ^ 13 
-gCs - gvr Ca + ^C3 - 



270 

17 



7r" + 



2137 



9 



-7r^ + 



24131 



72 



1360 170, 4 508 2, 2^ 1300 2, 787^ 
+ + — log^ 2 + — TT^ log' 2 - -^tt' log 2 - —(5 



340 2> 23311 . 20429 4 



+ 



1024 128, 4^ 256 3, 2. 
■04 + — log 2 + — tt' log' 2 - 



8705 
TO8 
1504 



1656817 



-TT 



1944 
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9 



9 

CpTpni 



tt' log 2 



+ 



1096, 916 4 904 , 1120 

H C3 H H 

3 135 9 9 

76804 + 32 log^ 2 - 327r' log' 2 + I^tt' log 2 



+ 



466 



+ 



9 

512 



C3 + 



124 
45" 



-TT 



8848 
128 



-TT 



9 

16811 
54 



752 



64 . _ 

04 log^ 2 - — TT^ log^ 2 + — TT^ log 2 

3 9^ 9 ^ 9 ^ 



280^ 152 4 52 2 111791 
"9"^^ + T35^ +^43" 



CaTfUi 
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,'8, 2461\ ,„^, , 64 . 64 2 5120 o, 

+ I 3C3 - ^ 1 «^ (/^) - 512a4 - y log^ 2 + y TT^ log2 2 + -^tt^ log 2 

, 44 2> 2837, 136 . 36268 2 1006271 ^ 

- 6OC5 + — vrYs Cs vr^ + C^Tp 

g S3 9 ^•^ 45 81 81 J 

/224 304 , 11534\ 2 2 ^ /208 , 18884\ 2 
V 9 ^ 27 243 y ^ ' V 9 81 y ^ 

/448 128 2 16850\ 2 

(here 04 = Li4(l/2)). Gauge dependence first appears at three loops, as in 
Zq [8]. The requirement of finiteness of the renormalized matching coeffi- 
cient ([S]) at e — )■ has allowed the authors of [S] to extract Zq from their 
result for Zq. 

It would not be too difficult to take into account a lighter massive flavour, say, 
rric 7^ in 6-quark HQET. Zq is no longer equal to 1, but is known at three 
loops [H]; Zq contains two scales, and is a non-trivial function of mc/mb [T5] . 
Both Zq and Zq have no smooth limit at rric — > 0, but the discontinuity 
cancels in the ratio (jl]). 

Now let's consider z{fi) in the large-/3o limit (see Chapter 8 in [16] for a ped- 
agogical introduction): 



.0 p V 2/* aft; ■ PoJo ' ' Vftv ' 

where (3 = /3oas/(47r), 7 = 7oa;s/(47r) + ■■■ (differences of n^-flavour and 
{rii + l)-fiavour quantities can be neglected at the l/(3o order). The difference 
of the QCD and HQET anomalous dimensions 'J = 'Jq — and the Borel 
image S{u) can be expressed as 

,(,)._4^,_,.0). 5M^ ^'°-">-^'»-°' . (9) 

Po u 

where the function F{e,u) has been calculated in [I] (see also [TB]): 

r(3-M-£) ^ ^ 

X (3-2£)(1-m)(1 + m-£), (10) 
D{e) = 6e^^"r(l + e)B{2 - e,2 - e) = 1 + + ■ ■ ■ 

The anomalous dimension difference [1] is gauge invariant at this order: 

.(B) = 2r ^ {i + m + m . .... 

^/3o5(2 + /3,2 + /3)r(3 + /3)r(l-/3) ' ^ ' 



Fie,u) = -2CFi^ 
\mj 
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the Borel image is [T7fT6 



S{u) 



JL+5/3} 



„ r(.x)r(i-2^) 



r(3 



— u 



u 



1 

2m 



(12) 



This Borel image has infrared renormalon poles at each positive half-integer 
u and at M = 2. Therefore, the integral in (|8]) is not well defined. Comparing 
its residue at the leading pole -u = 1/2 with the residue of the static-quark 
self-energy at its ultraviolet pole -u = 1/2 [18j, we can express the renormalon 
ambiguity of z{jj) as 



3AA 
2 m 



(13) 



(A is the ground-state meson residual energy). This ambiguity is compensated 
in physical matrix elements by ultraviolet renormalon ambiguities in the lead- 
ing 1/m correction (matrix elements of both local and bilocal dimension-5/2 
operators), see [17]. 

The matching coefficient is gauge invariant at the order l//3o- Expanding 7(/3) 
and S{u) and integrating, we obtain 



^(^) = l-C^^iM^3L + 4 + 



+ 
+ 



47r 



3 o 19, 2 113 

-L^ H L + TT^ H 

2 2 8 



/3o«s(/i) 



Atx 



167\ , ^ 19 2 5767 

L + 14(3 H H 

6 / ^3 216 



'/3oas(/i) 



An 



.4 19, s ( 2 167\,2 / > 2 2903 \ , 

+ -L' + f 37r2 + _ j + (36(3 + I^t^^ + ^ j ^ 



71 . 467 _ 167 2 



+ — vr^ + Cs + — 

^ 40 4 6 



103933 
1728 



47r 



+ 



(14) 



Thus we have confirmed the contributions with the highest power of ni in each 
term in ([7]), and predicted such a contribution at a^. 

Numerically, the matching coefficient ([7]) in the Landau gauge at fi = m and 
rii = 4 can be written as 



, , 4aW(m) , ^ , (at\my 

z(m) = 1 - - ' ^ ' - (1.9996/3o - 4.5421) ' ^ ' 



3 TT 



TT 



(2.2091/3o' + 5.1153/3o - 61.5397) |^^^5!!M 



- (3.3755/3^ + 



'a^^^(m) 



+ 



TT 



6 



^ 4 a(^)(m) 
3 



(16.6629-4.5421) 



- (153.4076 + 42.6271 



- (1953.4013 + ■ 



61.5397) 




Aaf\m) 
~ 3 



- 12.1208 




(15) 



- 134.4950 



TT 



'1953.4013 + 



(/So is for n; = 4 flavours). Naive nonabelianization [T] works rather well at 
two and three loops (in the latter case the C(/3o) and 0{1) terms partially 
compensate each other, similarly to [2]). Therefore, we can expect that the 
estimate of the term is also reasonably good. Numerical convergence of the 
series is very poor; this is related to the infrared renormalon at m = 1/2. 

Now let us consider the relation between the MS renormalized electron fleld 
in QED and the Bloch-Nordsieck electron fleld ^19j. The bare matching co- 
efficient Zq = Z^p is gauge invariant to all orders [20f2T|8] . In the Bloch- 
Nordsieck model, due to exponentiation [22], \ogZ^ = (3 — a(0))a(0)/(47r£) 
(where the 0-flavour a^^'' is equal to the on-shell a ~ 1/137). In the full QED, 
\ogZ^ = —a^^'^a^^y {Ane) + (gauge-invariant higher terms), see the Appendix. 
The gauge dependence cancels in \og{Z^/Z^) because of the QED decoupling 



relation a^^^a^^'' 



. Therefore, the renormalized matching coefficient 



z{iJ,) in QED is gauge invariant to all orders. We obtain 



z{fi) =1 - (3L + 4) 



a 



'9,2 31 2, o 55 2 5957\/a\2 

2^^ + yL + 16.^1og2-24C3-y.^ + ^ U 



2 19363\ ^ 



128, 4^ „ 2o 11 ''92 2 



+ 102404 + — log* 2 - 647r^ log^ 2 



3 

. ^ 9494 



72 J 
7r'log2 + 20C5 



104 

Is" 



259133 2 249887 



405 



324 



a 



V 



AttJ 



+ 



(16) 



In conclusion: we have derived the QCD/HQET matching coefficient for the 
heavy-quark fleld with three- loop accuracy ([7]), and the all-orders result in the 
large-/3o limit. The corresponding QED coefficient (fT6ll is gauge invariant. 
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Appendix. Gauge dependence of and 7^ in QED. 



The electron propagator S{x) is related to the Landau-gauge propagator Sl{x) 
(ao = 0) by the formula [20] 



S{x) = ^i(x)e-'^o(A(-)-A(o)) ^ (17) 

where A(x) is the Fourier image of A(/c) = ao/(A;^)^, and A(0) = in dimen- 
sional regularization. The electron field renormalization does not depend on its 
mass. For simplicity, we shall consider the massless electron, whose propagator 
has a single Dirac structure: 

S{x) = 5o(x)e'^(^) , (18) 
where Sq{x) is the li- dimensional free massless electron propagator. Then 

a{x) = a,{x) + ao^^ r(-.) , (19) 

where the Landau-gauge ai starts from Cg. Re-expressing a via renormalized 
quantities we have 

a{x) = a,{x) + a(/i)^ [^^) e^"^r(-.) . (20) 

This must be equal to \ogZ^ + a^, where log Z^(a(yu), a(/i)) contains only 
negative powers of e, and the renormalized — only non-negative. Therefore, 

a 

log ZJ a, a) = log Ziia) - a- — . (21) 

In QED (ilog(a(yu)a(/i))/(ilog/i = —2e exactly, anc{3 

a 

7^(a,a) = 2a— + 7L(a) , (22) 
47r 

where the Landau-gauge 7L(a) starts from a^. 



^ I was informed by D.J. Broadhurst and D.V. Shirkov that this result has been 
proved in some Russian article in the second half of 50s. I am grateful to them for 
discussing this question; unfortunately, I was unable to find this article. 
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